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Abstract 

The effect of nearest-neighbor coupling on the thermodynamic and dynamical properties of the 
ferromagnetic Hamiltonian Mean Field model (HMF) is studied. For a range of antiferromagnetic 
nearest-neighbor coupling, a canonical first order transition is observed, and the canonical and 
microcanonical ensembles are non-equivalent. In studying the relaxation time of non-equilibrium 
states it is found that as in the HMF model, a class of non-magnetic states is quasi-stationary, with 
an algebraic divergence of their lifetime with the number of degrees of freedom N. The lifetime of 
metastable states is found to increase exponentially with TV as expected. 
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1. INTRODUCTION 



, this result has been confirmed 
, which pointed 



Hamiltonian systems with long-range interactions exhibit several peculiarities, both in 
their equilibrium properties and in the dynamical behavior [lj]. At equilibrium, it is often 
found that different statistical ensembles can be non-equivalent even in the thermodynamic 
limit. First discussed theoretically by Hertel and Thirrine 
in analytical and numerical studies (see, e.g., Refs. 

m 

out that ensemble inequivalence is usually associated with the presence of a first order phase 
transition in the canonical ensemble. 

Attention has also been focussed on the dynamics, in particular on the relaxation to equi- 
librium of a system initially prepared in a state which is far from equilibrium. Many studies 
have been devoted to the Hamiltonian Mean Field (HMF) model (l2j, a system of classical 
XY-rotators with infinite range ferromagnetic couplings. This model exhibits a second order 
transition to a magnetized state. It has been shown numerically in microcanonical simu- 
lations that for a large class of non-equilibrium initial conditions the model exhibits fast 
relaxation to another non-equilibrium state, in a fashion similar to the violent relaxation 



of astrophysical systems 
equilibrium takes p 
a power law [l2j, |l5 
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14j . After this transient the relaxation to Boltzmann-Gibbs 



ace in times that diverge with the number N of rotators, according to 
The stability of these long-lived states, that have been called 



"quasi-stationary" states, has recently been analyzed in terms of Vlasov and Fokker-Planck 
equations [ill I2H 2Jj • It should be emphasized that these quasi-stationary states are 
not metastable states in the thermodynamical sense, since they do not correspond to local 
maxima of the entropy (in the microcanonical ensemble) or local minima of the free energy 
(in the canonical ensemble). A similar effect has recently been observed in the Ising model 
with long and short range interactions, where the lifetime of the quasi- stationary states 
diverges logarithmically with N (2^]. 

In this work we consider the robustness of the quasi- stationary states with respect to 
perturbations of the long-range Hamiltonian. This problem is of interest if one wants to 
ascertain the relevance of these states in realistic long-range systems, where a purely long- 
range Hamiltonian could be perturbed by external forces or by the presence of additional 
short-range terms in the potential energy. In particular, we generalize the HMF model 
to include nearest-neighbor interaction between the rotators, in addition to the long range 
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interaction. 

The paper is organized as follows: In Section |21 we introduce our model; in Section |3] we 
summarize the equilibrium properties of the model as obtained within the canonical and 
the microcanonical ensembles. In Sections 0] and we report our numerical results on the 
relaxation time of quasi-stationary and metastable states, respectively. Concluding remarks 
are given in Section |U1 



2. THE MODEL 



We study a system of N rotators characterized by phase variables 0j, whose Hamiltonian 
is given by: 

N 2 j N N 

The N rotators are placed on a one-dimensional lattice with periodic boundary conditions 
(#jv+i = 6\). The choice of a one-dimensional lattice is made for mathematical convenience, 
as it will be explained in the following Section. The parameter K is the short-range (ac- 
tually nearest-neighbor) coupling between rotators. For K = this system reduces to the 
HMF model, that, for positive J, has a second order ferromagnetic transition at the critical 
temperature T c = 0.5J, corresponding to the critical energy density e c = 0.75J 
in this case the topology of the lattice is irrelevant. The transition is characterized by the 
value of the order parameter given by the magnetization m: 



12]. 



Clearly 



\ 



N \ * / N 



^TcosM + PT sin^ =Jml + ml. (2) 



,1=1 / \i=l 

In the iV — > oo limit, m is positive below the critical energy, approaching zero for e — > e c ; 
while m is identically zero above e c . It has been shown that for the HMF model microcanon- 
ical and canonical ensembles are equivalent ^j] ; this equivalence had been discussed before 
for a general class of magnetic systems with long range interactions j^j. Microcanonical 
simulations in which the system is initially prepared in an homogeneous state, i.e. with the 
angles Qi uniformly distributed between and 27r (and thus m w at time t — 0) and the 
Pi's uniformly distributed between — p and po (with p determined by the energy) show the 
following: if the energy per particle is in a certain range below the critical energy, then the 



system remains for a long time in this non-magnetic state , an d it approaches the equilibrium 
magnetized state only after a time that diverges like iV L7 In the same subcritical energy 
range it is found that if the 9^s are initially all equal, with the p^s uniformly distributed, 
there is a fast relaxation to a non-magnetic state, and again the equilibrium magnetized 
state is reached after a time diverging with N, although with a smaller power law 17]. 

In the following sections we study the equilibrium properties of system (0) for general 
values of K, and consider the dynamical properties at small values of K, when the system 
can be thought of perturbed HMF. 



3. EQUILIBRIUM PROPERTIES 



The canonical partition function and the microcanonical entropy of the model can be 
obtained by a straightforward generalization of the corresponding calculations of the HMF 



model 
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121 ] . We will not report here the full details of these calculations, which we defer 



to a longer paper, but rather restrict ourselves to a brief sketch of the derivation. Without 
loss of generality we take J = 1 in ((TJ. 

We begin by considering the free energy per particle f(/3), obtained from the calculation 
of the canonical partition function. Applying the Hubbard- Stratonovich transformation, 
f((3) can be written as 



max 

m 



(3(1 + m 2 ) 2tt 

+ In A(/3m, (3K) + In — 

2 p ^ 



(3) 



In this expression X(z, a) is the largest eigenvalue of the transfer matrix given by the sym- 
metric integral operator T defined as 



(T0)(0) = / dO'exp 



-z(cos 9 + cos 9') + a cos(6> - 9') 



(4) 



The energy density e{(3) is computed from Q using standard expressions. The microcanon- 
ical entropy as a function of e can be computed from the argument of maximization in 
(J2J) using the recipe introduced in Ref. |7j, performing first a minimization over the inverse 
temperature (3 before maximizing over the order parameter m 

2vr (3(1 + m 2 ) 



s e 



max mm 

m p 



+ In \(p m , (3K) 



= max s(e, m) 



(5) 



The value of (3 corresponding to the solution of the optimization problem (jSJ) gives the 
microcanonical temperature as a function of e. The phase-diagram in both the canonical 
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FIG. 1: The canonical and micro canonical (K,T) phase diagram. In the canonical ensemble the 
transition is continuous (bold solid line) down to the tricritical point CTP where it becomes first 
order (dashed line). In the microcanonical ensemble the continuous transition coincides with the 
canonical one at large K (bold solid line). It persists at lower K (light solid line) down to the 
tricritical point MTP where it turns first order, with a branching of the transition line (dotted 
lines). 

and the microcanical ensembles can be directly derived from the properties of free energy and 
entropy. For positive (i.e., ferromagnetic) nearest-neighbor couplings the system behaves in 
a way which is analogous to the HMF (K = 0) case: there is a second order ferromagnetic 
transition (at a K dependent critical temperature and energy), ensembles are equivalent and 
the temperature-energy relation (so-called caloric curve) is qualitatively similar to that of 
the HMF model. It is found that there exists a range of negative values of the short-range 
coupling K for which the system exhibits a first order phase transition in the canonical 
ensemble. This is a result of the competing effect of the mean-field potential and the 
antiferromagnetic short-range coupling. As expected, the canonical and the microcanonical 
ensembles are non-equivalent in this range of coupling. The (K, T) phase diagram of the 
model is given in Fig. ^ For clarity the plot is restricted to the interesting negative K 
region. Starting from positive values and decreasing K, the canonical critical line ends at 
the tricritical point (CTP) K = K\ ~ —0.168, T = T\ ~ 0.273, where the canonical phase 
transition becomes first order. The dashed line in Fig. Q shows the canonical transition 
temperature for K < K\. This line ends at K — —0.25, T = 0. For smaller values of K 
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FIG. 2: Two examples of s(e,m): (a) Local entropy minimum for a m = quasi-stationary state; 
(b) Local entropy maximum for a metastable m = state. 

the system is disordered at any temperature and there is no transition. This can be easily 
understood by computing the energies of the staggered antiferromagnetic m = state and 
that of the fully magnetized m = 1 state. For K < —1/4 the former is always favoured, 
even at zero temperature. In the microcanonical ensemble the critical line extends down to 
the microcanonical tricritical point (MTP) located at K = K2 ~ —0.182, T = T2 ~ 0.234, 
where the transition becomes first order. Below this point, a temperature jump between the 
two microcanonical temperatures given by the dotted lines appears. Again, for K < —0.25 
there is no transition. The two ensembles are therefore non-equivalent for —0.25 < K < K\. 



4. QUASI-STATIONARY STATES 

The equilibrium magnetization at a given energy e is obtained by maximizing the entropy 
s(e, m) with respect to m (see Eq. (|5|)). Local maxima correspond to metastable states, that 
we consider in the next Section. For small K values (both positive and negative) the entropy 
dependence on magnetization shows only a global maximum, at the equilibrium value of m. 
A typical entropy curve is given by the one shown in Fig. |2Ja), which has K — (HMF 
model); the entropy has a local minimum at m = 0. We are primarily interested in this case 
since we want to study the difference in the dynamical behavior induced by a perturbation 
of the HMF system, for which the long lifetime of quasi- stationary states, with m close to 
0, is strictly of dynamical origin Q]. We consider here only the case of initial homogeneous 



6 





5000 10000 15000 20000 
t 

FIG. 3: Evolution of the magnetization for quasi-stationary states with increasing number of 
rotators. Here K = 0.05 and e = 0.7115. From left to right N = 1024 (100 runs), N = 2048 (30 
runs), N = 4096 (10 runs), N = 8192 (10 runs), N = 16384 (2 runs), N = 32768 (2 runs). 

non-magnetic states, with 9 and p both uniformly distributed. 

We performed runs at different K values. For each K we have chosen a value of the 
energy density e below the critical energy, e = e c (K) — Ae, keeping Ae = 0.06 fixed. This 
distance is the same as for the simulations typically performed at K = (HMF), where 
e = 0.69 is studied, with the critical energy density being at 0.75 We have performed 
runs for a number of rotators iV ranging from 1024 to 32768. 

It turns out that the escape from non-magnetic quasi- stationary states takes place after a 
time that fluctuates from run to run (corresponding to different realizations of the uniform 
initial distributions of 9 and p), although these fluctuations are much smaller than those of 
the case of the thermodynamic metastable case discussed in the following Section. In Fig. |3] 
we plot the instantaneous magnetization as a function of time for the case K = 0.05 and 
e = 0.7115, and for the increasing N values from left to right. These plots are obtained 
by averaging over several runs (the number of these runs for each N value is reported in 
the caption). It is worth mentioning that for a given N the lifetime of the quasi-stationary 



. 3). 



To define 



states is considerably shorter than that of the HMF model (see, e.g., Ref. 
this lifetime we have chosen the time at which, during the rise of the magnetization towards 
the equilibrium value, the magnetization has a value equal to the average between the 
plateau value in the quasi-stationary state and the equilibrium value. In order to smooth 
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FIG. 4: Growth of the lifetime of quasi-stationary states with iV for different values of K. From 
top to bottom the data refer to K = O^K = 0.0025, K = 0.05, K = 0.1. Data for K = 0, with r 
diverging as TV 1 ' 7 , are taken from Ref. 18]. 



out fluctuations, we used a running average of the magnetization over a sliding time window 
of a conveniently chosen size. This definition of the lifetime is practically equivalent to that 
used in Ref. [18] through a more sophisticated fitting procedure. In Fig.HJwe plot the lifetime 
r as a function of N, in a log-log plot for different K values. It is evident that the lifetime 
r increases with N, possibly diverging in the limit iV — > oo. This suggests that quasi- 
stationarity persists when short range interactions are added to the HMF model. A naive 
fit of the curves to the form r ~ iV 7 with i^-dependent power- law exponents 7 (if) yields 
7(0.0025) ~ 0.65, 7(0.05) ~ 0.56, 7(0.1) w 0.45. Thus the observed exponent 7 decreases 
with increasing K. For K > 0.1 we do not observe a plateau region in the magnetization at 
m K and one cannot reasonably speak of quasi-stationarity. Clearly one cannot rule out 
the possibility that the seemingly K dependence of 7 is just a finite size effect, which are 
evidently present in the up-bending of the smaller K value curve in Fig. 0] and the down- 
bending of the other two. More extensive numerical simulations are necessary in order to 
assess the law of divergence with N. 
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FIG. 5: Evolution of the magnetization for the metastable situation in which the system with e 
and K as in Fig. 2(b) is started with m = 0. 

5. METASTABLE STATES 



In our system, thermodynamic metastable states are found for negative K values, for 
which the microcanonical transition is first order. In Fig. Efb) we show an example of an 
entropy function with a local maximum at m — 0, while the global maximum is at m ^ 0. 
The evolution of this m = state was studied. In Fig. we show a typical plot of the 
istantaneous magnetization m as a function of time, for a system with iV = 150 rotators. 
After a transient time spent in the metastable state, there is a sudden transition to the 
equilibrium magnetization. In Fig. El we plot the lifetime r of the metastable state as a 
function of N. It should be pointed out that an average over many runs is necessary in 
order to obtain a sufficiently accurate determination of this value, since the differences from 
run to run are not small. We find that r ~ e NAs , where As is the "entropy barrier", i.e. the 
difference between the entropy at m = and the one at the minimum (see Fig. Efb)). Thus 
the transition is determined by the fluctuations at finite N, and the exponential divergence 
of the lifetime is governed by the entropy barrier. Analogous results have been found for 



other systems 



241- 
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FIG. 6: Exponential growth of the lifetime of the metastable state of Fig. 3 with the number of 
rotators. 

6. DISCUSSION 



The model studied in this work offers another example of a long-range system with 
inequivalence of ensembles. When the long-range and the short-range interactions have a 
competing character (i.e., when K is negative), we found that the microcanonical and the 
canonical ensembles yield different results, in a range of K values. 

The existence of quasi-stationary states seems to be a peculiar feature of long-range 
systems, and in this work our aim was to study their robustness with respect to perturbations 
of the long-range nature of the interaction. We found that if the perturbation is sufficiently 
small the divergence with N of the lifetime of these states is preserved. This fact was not 
obvious a priori, since one could have argued that even small perturbations could have given 
rise to finite lifetimes also in the N —* oo limit. One could have hypothesized the finiteness 
of the lifetime on the basis of the fact that, when the rotators have also a short-range 
coupling, the 1/N expansion of the Vlasov equation, suitable for the HMF case j^], is no 
longer possible. It has to be inferred that, in spite of this fact, an approximate description 
of the system with a kinetic equation possessing out of equilibrium stable stationary states 
would still be possible. Its determination is one of the possible directions of development of 
this work. It is worth noting that in the case of the Ising model with both long and short 
range interactions, quasistationarity was found to exist, but with a universal logarithmic 
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divergence of the lifetime with iV 22]. 



We have considered here only one class of initial conditions. As was found in the HMF 
model, other initial states may behave differently. It would be interesting to explore the 
dynamical behavior of other initial conditions. The model we have analyzed lives on a ID 
lattice, but we believe that some basic dynamical features observed here should extend to 
higher dimensions. 

It would of course be interesting to determine the range of validity, across different models, 
of the results presented in this paper; this could be relevant for the study of the dynamical 
properties, and in particular for the approach to equilibrium of open long-range systems, i.e., 
systems in interaction with external degrees of freedom. One could suppose that the open 
character of the system can be mimicked by a perturbation of the long-range Hamiltonian. 



Preliminary investigations of a simple magnetic system with these characteristics 20J give 
results with features that are similar to those presented in this work. 
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